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Abstract 

Data from e"'"e~ annihilation into hadrons collected by the JADE experiment at 
centre-of-mass energies between 14 GeV and 44 GeV are used to study the four- 
jet event production rate as a function of the Durham jet algorithm's resolution 
parameter ycut- The four-jet rate is compared to QCD next-to-leading order calcu- 
lations including resummation of large logarithms in the next-to-leading logarithmic 
approximation. The strong coupling measured from the four-jet rate is 

as(Mzo) = 0.1159 ±0.0004(stat.) ± 0.0012(exp.) ± 0.0024(had.) ± 0.0007(theo.) , 

in agreement with the world average. 
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1 Introduction 



The annihilation of electrons and positrons into hadrons allows precise tests of Quantum 
Chromodynamics (QCD). Many observables have been devised which provide a conve- 
nient way of characterizing the main features of such events. Multijet event rates are 
predicted in perturbation theory as functions of the jet-resolution parameter, with one 
free parameter, the strong coupling a^. Events with four quarks in the final state, qqqq, 
or two quarks and two gluons, qqgg, may lead to events with four-jet structure. In leading 
order perturbation theory, the rate of four-jet events in e^e~ annihilation is predicted to 
be proportional to a^. The strong coupling can be measured by determining the four-jet 
event production rate and fitting the theoretical prediction to the data. 
Calculations beyond leading order are made possible by theoretical progress achieved dur- 
ing the last few years. For multi-jet rates as well as numerous event-shape distributions 
with perturbative expansions starting at 0{as), matched next-to-leading order (NLO) 
and next-to-leading logarithmic approximations (NLLA) provide a satisfactory descrip- 
tion of the data over large kinematically allowed regions at many centre-of-mass energies 
[2-5]. 

First evidence for four-jet structure has been reported earlier by the JADE collabora- 
tion [6]. In addition multi-jet event production rates were measured and the three-jet 
rate was used to determine the the strong couphng as [7-9]. The ALEPH, DELPHI and 
OPAL collaborations published measurements of as based on the four-jet rate in the en- 
ergy range between 91 and 209 GeV [10-12]. The same theoretical predictions as used 
here were employed to determine the strong coupling as- 

In this analysis we use data collected by the JADE experiment in the years 1979 to 1986 
at the PETRA e''"e~ collider at DESY at six centre-of-mass energies spanning the range 
of 14-44 GeV. 



2 Observable 

Jet algorithms are applied to cluster the large number of particles of a hadronic event into 
a small number of jets, refiecting the parton structure of the event. For this analysis we 
use the Durham scheme [2]. Defining each particle initially to be a proto-jet, a resolution 
variable is calculated for each pair of proto-jets i and j: 

2mm{Ef, E]) 

Vij = ^2 ^(1 - cos %), (1) 

where Ei and Ej are the energies of jets i and j, cos 6ij is the cosine of the angle between 
them and E^[s is the sum of the energies of the detected particles in the event (or the 
partons in a theoretical calculation). If the smallest value of yij is less than a predefined 
value ^cut; the pair is replaced by a new proto-jet with four- momentum — Pi +Pj, and 
the clustering starts again. Clustering ends when the smallest value of Uij is larger than 
Ucnt, and the remaining proto-jcts are accepted and counted as finally selected jets. 
The next-to-leading order QCD calculation predicts the four-jet rate i?4, which is the 
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fraction of four-jet events as a function of as- It can be written in the following way [5]: 

R^iycnt) = = v'B,{y^t) + v'[C4ycut) + 3/2(/3o logx^ - 1) B^y^^t)] (2) 

where (T4-jct(ycut) is the cross-section for the production of hadronic events with four jets at 
fixed i/cut, CTtot the total hadronic cross-section, r] = as Cf/271 with Cp = 4/3 the colour 
factor from the SU(3) symmetry group, = [xj^fs with \i being the renormalization 
scale, ^/s the centre-of-mass energy, and /3o = (11 — 2n//3) with n/ the number of active 
flavours^. The coefficients i?4 and C4 are obtained by integrating the massless matrix 
elements for e+e^ annihilation into four- and five-parton final states, calculated by the 
program DEBRECEN 2.0 [5]. 

Eq. (2) is used to predict the four-jet rate as a function of ^cut- The fixed-order perturba- 
tive prediction is not reliable for small values of ^cut, due to terms of the form ccg ln"^(ycut) 
which enhance the higher order corrections. An all-order resummation of such terms in 
the NLLA is possible for the Durham clustering algorithm [2]. The NLLA calculation is 
combined with the NLO-prediction using the so-called modified R-matching scheme [13]. 
In the modified R-matching scheme the terms proportional to rf' and rf are removed from 
the NLLA prediction i^^LLA ^j^^ difference is then added to the NLO calculation F!^^^ 
with the result: 

^R-match _ ^NLLA ^ [^2^^^ _ ^NLLA^ ^ ^3^^^ _ ^^NLLA _ 3/3(5^ _ 5^^^^))], (3) 

where iJ^LLA ^^^^ (^nlla ^j^g coefficients of the expansion of i^^LLA ^ ^^2^ 
the term and 7/cut dependence have been suppressed for clarity. 



3 Analysis Procedure 
3.1 The JADE Detector 

A detailed description of the JADE detector can be found in [1]. This analysis relies 
mainly on the reconstruction of charged particle trajectories and on the measurement of 
energy deposited in the electromagnetic calorimeter. Tracking of charged particles was 
performed with the central tracking detector, which was positioned in a solenoidal magnet 
providing an axial magnetic field of 0.48 T. The central detector contained a large volume 
jet chamber. Later a vertex chamber close to the interaction point and surrounding z- 
chambers to measure the z-coordinate ^ were added. Most of the tracking information is 
obtained from the jet chamber, which provides up to 48 measured space points per track, 
and good tracking efficiency in the region |cos^| < 0.97. Electromagnetic showers are 
measured by the lead glass calorimeter surrounding the magnet coil, separated into a barrel 

^In this analysis the number of active flavours is set to five. 

•^In the JADE right-handed coordinate system the -\-x axis points towards the centre of the PETRA 
ring, the y axis pointed upwards and the z axis points in the direction of the positron beam. The polar 
angle Q and the azimuthal angle ^ are defined with respect to z and x, respectively, while r is the distance 
from the z-axis. 
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(jcos^l < 0.839) and two end-cap (0.86 < |cos^| < 0.97) sections. The electromagnetic 
calorimeter consisted of 2520 lead glass blocks with a depth of 12.5 radiation lengths in 
the barrel (since 1983 increased to 15.7 in the middle 20% of the barrel) and 192 lead 
glass blocks with 9.6 radiation lengths in the end-caps. 

3.2 Data Samples 

The data used in this analysis were collected by JADE between 1979 and 1986 and corre- 
spond to a total integrated luminosity of 195 pb^^. The breakdown of the data samples, 
average centre-of-mass energy, energy range, data taking period, collected integrated lu- 
minosities and the size of the data samples after selection of hadronic events are given in 
table 1. The data samples are chosen following previous analyses, e.g. [1,9,14-16] ^. 



average 
energy in GeV 


energy 
range in GeV 


year 


luminosity 
(Pb-^) 


selected 
events 


14.0 


13.0-15.0 


1981 


1.46 


1783 


22.0 


21.0 23.0 


1981 


2.11 


1103 


34.6 
35.0 


33.8-36.0 
34.0-36.0 


1981-1982 
1986 


61.7 
92.3 


14313 
20876 


38.3 


37.3-39.3 


1985 


8.28 


1585 


43.8 


43.4-46.4 


1984-1985 


28.8 


4374 



Table 1: The average centre-of-mass energy, the energy range, the year of data taking 
and the integrated luminosity for each data sample, together with the numbers of selected 
data events. 



3.3 Monte Carlo Samples 

Samples of Monte Carlo simulated events are used to correct the data for experimental 
acceptance, resolution and backgrounds. The process e"'"e~ hadrons is simulated using 
PYTHIA 5.7 [17]. Corresponding samples using HERWIC 5.9 [18] are used for systematic 
checks. The Monte Carlo samples generated at each energy point are processed through a 
full simulation of the JADE detector [19-22] , and reconstructed in essentially the same way 
as the data. In addition, for comparisons with the corrected data, and when correcting 
for the effects of hadronization, large samples of Monte Carlo events without detector 
simulation are employed, using the parton shower models PYTHIA 6.158, HERWIC 6.2 
and ARIADNE 4.11 [23]. Each of these hadronization models contains a number of 
tunable parameters; these were adjusted by tuning to previously published OPAL data 
at ~ 91 GeV as summarized in [24] for PYTHIA/JETSET and in [25] for HERWIG 
and ARIADNE. The data taken by the JADE experiment are well described by the Monte 

^The data are available from two versions of the reconstruction software from 9/87 and from 5/88. We 
use the set from 5 /88 as the default version and consider differences between the results as an experimental 
systematic uncertainty. 
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Carlo simulated events [16], with only the centre-of-mass energy changed and set to the 
respective value of the event. 

3.4 Selection of Events 

The selection of events for this analysis aims to identify hadronic event candidates and 
to reject events with a large amount of energy emitted by initial state photon radiation 
(ISR). The selection of hadronic events is based on cuts on event multiplicity to remove 
leptonic final states and on visible energy and longitudinal momentum imbalance to re- 
move radiative and two-photon events, e"'"e~ — > e"'"e'"-|- hadrons. The cuts are documented 
in [8,26,27] and summarized in a previous publication [14]: 

Standard criteria are used to select good tracks and clusters of energy deposits in the 
calorimeter for subsequent analysis. For charged particle tracks the pion mass was as- 
signed while for clusters zero mass particles are assumed. Charged particle tracks are 
required to have at least 20 hits in r — and at least 12 hits in r — 2; in the jet chamber. 
The total momentum is required to be at least 100 MeV. Furthermore, the point of clos- 
est approach of the track to the collision axis is required to be less than 5 cm from the 
nominal collision point in the x — y plane and less than 35 cm in the 2;— direction. 
In order to mitigate the effects of double counting of energy from tracks and calorimeter 
clusters a standard algorithm is adopted which associates charged particles with calorime- 
ter clusters, and subtracts the estimated contribution of the charged particles from the 
cluster energy. Clusters in the electromagnetic calorimeter are required to have an energy 
exceeding 0.15 GcV after the subtraction of the expected energy deposit of associated 
tracks. From all accepted tracks and clusters the visible energy E^is = J2iEi, momentum 
balance Pbai = I J2iPz,i\/ E^is and missing momentum Pmiss = I SiPil are calculated. 
Hadronic event candidates are required to pass the following selection criteria: 

• The total energy deposited in the electromagnetic calorimeter has to exceed 1.2 GeV 
(0.2 GeV) for < 16 GeV, 2.0 GeV (0.4 GeV) for 16 < < 24 GeV and 3.0 GeV 
(0.4 GeV) for ^/s > 24 GeV in the barrel (each endcap) of the detector. 

• Events with exactly four tracks with configurations where three tracks are in one 
hemisphere and one track is in the opposite hemisphere are rejected. 

• At least three tracks have to have more than 24 hits in r — and a momentum in 
the X — y plane larger than 500 MeV; these tracks are called long tracks. 

• The visible energy has to fulfill £'vis/-\/s > 0.5. 

• The momentum balance has to fulfill pbai < 0.4. 

• The missing momentum has to fulfill Pmiss/ a/s < 0.3. 

• The z-coordinate of the reconstructed event vertex has to lie within 15 cm of the 
interaction point. 

• The polar angle of the event thrust axis [28] is required to satisfy | cos(^^t)| < 0.8 in 
order that the events be well contained within the detector acceptance. 
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• The number of good charged particle tracks is required to be greater than three 
reducing t+t~ and two-photon backgrounds to a neghgible level. 

The numbers of selected events for each centre-of-mass energy are shown in table 1. 
3.5 Corrections to the data 

All selected tracks and electromagnetic calorimeter clusters are used in the evaluation of 
the four-jet rate. The four-jet rate distribution as a function of the jet resolution ycut 
after all selection cuts applied is called the detector-level distribution. 
In this analysis events from the process e"^e~ bb systematically bias our results, 
since especially at low centre-of-mass energies the large mass of the b quarks and of 
the subsequently produced and decaying B hadrons will influence the four-jet rate dis- 
tribution [10-12]. The QCD predictions are calculated for massless quarks and thus we 
choose to correct our data for the presence of bb events. About 1/11 of all qq events 
are bb events. The expected number of bb background events rji is subtracted from the 
observed number of data events Ni at each t/cut bin ^- The effects of detector acceptance 
and resolution and of residual ISR are then accounted for by a multiplicative correction 
procedure. 

Two distributions are formed from Monte Carlo simulated signal events; the first, at the 
detector-level, treats the Monte Carlo events identically to the data, while the second, at 
the hadron-lcvel, is computed using the true four- momenta of the stable particles^, and 
is restricted to events where V^, the centre-of-mass energy reduced due to ISR, satisfies 
y/s — y/s' < 0.15 GeV. The ratio of the Monte Carlo distribution taken at the hadron-level 
to the distribution taken at the detector- level for each y^ut bin i, C'detector^ -g ^g^^^ ^^g ^ 
correction factor for the data. This finally yields the corrected number of four jet events 
s-t Z/cut bin i, = C'detector . j^jY^ — rji). The hadron-level distribution is then normalized at 
each i/cut bin i by calculating i?4_j = Ni/N, where the is the corrected number of events 
selected at hadron-level. We observe some disagreement between the detector corrections 
calculated using PYTHIA or HERWIG at low centre-of-mass energies while at larger ^/s 
the correction factors agree well within the regions chosen for comparison with the the- 
ory predictions. The difference in detector corrections is evaluated as an experimental 
systematic uncertainty. The numerical results of the four-jet rate at hadron-level at the 
different energy points are summarized in Tables 5 and 6. 



4 Systematic Uncertainties 

Several sources of possible systematic uncertainties are studied. Uncertainties originating 
from massless quark calculations arc not considered, since contributions to the four-jet 
rate from B hadrons are subtracted at detector-lcvcl. For each variation of parameters 
the difference of the resulting value of as with respect to the default value is taken as a 
systematic uncertainty. The default value of as is determined with the standard event 

^AU charged and neutral particles with a lifetime longer than 3 x 10~^° s are treated as stable. 
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selection and the correction procedure using PYTHIA. The systematic uncertainty is 
taken to be symmetric around the default value of as- 

4.1 Experimental uncertainties: 

Contributions to the experimental uncertainties arc estimated by repeating the analysis 
with varied cuts or procedures. For each systematic variation the value of as is determined 
and then compared to the result of the standard analysis (default value). 

1. In the standard analysis the reconstruction software from 5/88 is used. As a varia- 
tion a different reconstruction software from 9/87 is used. 

2. In the default method the estimated minimum ionizing energy from tracks associated 
with electromagnetic calorimeter clusters is subtracted from the cluster energies. As 
a variation all accepted tracks and all uncorrected electromagnetic clusters are used. 

3. The thrust axis is required to satisfy |cos(^t)| < 0.7. With this more stringent 
cut events are restricted to the barrel region of the detector, which provides better 
measurements of tracks and clusters compared to the endcap regions. 

4. Instead of using PYTHIA for the correction of detector effects as described in sec- 
tion 3.5, events generated with HERWIG are used. 

5. The requirement on missing momentum is dropped or tightened to Pmiss/v^ < 0.25. 
The larger deviation from the default value is taken as a systematic uncertainty. 

6. The requirement on the momentum balance is dropped or tightened to pbai < 0.3. 
The larger deviation from the default value is taken as a systematic uncertainty. 

7. The requirement on the number of long tracks is tightened to Mong > 4. 

8. The requirement on the visible energy is varied to E^^s/ > 0.45 and Ey^s/ y/s > 
0.55. The larger deviation from the default value is taken as a systematic uncer- 
tainty. 

9. The fit range is changed. Two different cases are considered. First the fit range is 
reduced by one data point at each edge of the standard fit range. Second the fit 
range is extended by one data point at each edge of the standard fit range. The 
larger deviation from the default fit is taken as a systematic uncertainty. In order 
to take statistical fluctuations into account, the deviation is calculated using the 
average deviation of a fit applied to 50 Monte Carlo samples. 

10. The amount of subtracted bb background is varied by ±5% of its nominal value of 
about 1/11 to cover uncertainties in the estimation of the background fraction in 
the data. The larger deviation from the default value is taken as the systematic 
uncertainty. 

All contributions listed above are added in quadrature and the result is quoted as the 
experimental systematic uncertainty. The dominating effects are the use of the different 
data versions and the different correction for detector effects. 
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4.2 Hadronization: 



The uncertainties associated with the hadronization correction (see section 5.2) are as- 
sessed by using HERWIG and ARIADNE instead of the default hadronization correction 
using PYTHIA. The larger change in resulting from these alternatives is taken to define 
the symmetric hadronization systematic uncertainty. 

4.3 Theoretical Uncertainties: 

The theoretical uncertainty, associated with missing higher order terms in the theoretical 
prediction, is assessed by varying the renormalization scale factor x^. The predictions of 
a complete QCD calculation would be independent of x^, but a finite-order calculation 
such as that used here retains some dependence on x^. The renormalization scale factor 
Xf^ is set to 0.5 and one. The larger deviation from the default value of as is taken as 
systematic uncertainty. 

5 Results 

5.1 Four- Jet Rate Distributions 

The four-jet rates for the six centre-of-mass energy points after subtraction of bb back- 
ground and correction for detector effects are shown in figures 1 and 2. Superimposed 
are the distributions predicted by the PYTHIA, HERWIG and ARIADNE Monte Garlo 
models. Towards large ijcut values (right to the maximum of the distribution) the decrease 
of the four-jet rate corresponds to the migration and classification to three- and two-jet 
events. Towards smaller T/cut values (left to the maximum of the distribution) the decrease 
corresponds to the migration and classification to five or more jet events, i.e. towards the 
higher order QGD and non-perturbative or hadronization region. In order to make a 
more clear comparison between data and models, the inserts in the upper right corner 
show the differences between data and each model, divided by the combined statistical 
and experimental error at that point. The sum of squares of these differences would, in 
the absence of point-to-point correlations, represent a between data and the model. 
However, since correlations are present, such values should be regarded only as a rough 
indication of the agreement between data and the models. The three models are seen to 
describe the data well. 

5.2 Determination of cxs 

Our measurement of the strong coupling as is based on fits of QCD predictions to 
the corrected four-jet rate distribution, i.e. the data shown in figures 1 and 2. The 
theoretical predictions of the four-jet rate using the combined 0(Q;g)+NLLA calculation 
as described in section 2 provide distributions at the parton-level. In the Monte Garlo 
simulation the parton-level distributions are obtained from the partons after the parton 
shower has stopped, just before the hadronization. In order to confront the theory with 
the hadron-level data, it is necessary to correct for hadronization effects. This is done by 
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calculating the distributions at both the hadron and the parton-level using PYTHIA and, 
as a cross-check, with the HERWIG and ARIADNE models. The theoretical prediction 
is then multiplied by the ratio of the hadron- and parton-level distributions. 

The differences between the models are considered as a systematic uncertainty. 
A x^-value at each energy point is calculated using the following formula: 

n 

= E(^4, - R{asyiniV{R,r%{R,j - Riasfin (4) 



theo 
4,i 



where the indices i and j denote the ycut bins in the chosen fit range and the R{as) 
are the predicted values of the four-jet rate. The four-jet rate as a function of 7/cut is an 
integrated distribution and therefore a single event can contribute to several ijcnt bins and 
neighboring 7/cut bins are correlated. The covariance matrix Vij is calculated using the 
matrix Wij determined from four-jet rate distributions calculated from Monte-Carlo at 
the hadron-level as follows: 



n n n 

.k=l fc=l fc=l 



(5) 



where is the average jet rate for a given y^ut bin i and sample k, and n is the number 
of subsamples. Subsamples are built by choosing 1000 events randomly out of the set 
of all generated Monte Carlo events. A single event may be included in several subsam- 
ples, but the impact on the final covariance matrix is expected to be very small and, 
therefore, is neglected [29]. For every centre-of-mass energy point 1000 subsamples are 
built. The matrix Wij is then used to determine the correlation matrix, pij = Wij/ai(fj, 
with (fj = \JWii. The covariance matrix Vij{R^ used in the fit is then determined 
using the statistical error Ui of the data sample at data point i and the correlation matrix 

The value is minimized with respect to as for each centre-of-mass energy point sepa- 
rately. The renormalization scale factor x^, as discussed in section 2, is set to one. 
The fit ranges are determined by requiring that the hadronization corrections be less than 
50% and the detector corrections be less than 50% in the fit range. In order to exclude 
the non-perturbative region we require • i/cut to be larger than 2 GeV. In the Durham 
scheme this value corresponds to the minimal transverse momentum of the pair of proto- 
jets with respect to each other. The fit range is 0.0209 < ycut < 0.0495 for data taken 
at 14 GeV, 0.0088 < y^ut < 0.0495 for data taken at 22 GeV, 0.0037 < y^nt < 0.0279 for 
data taken at 34.6 and 35 GeV, 0.0028 < y^nt < 0.0279 for data taken at 38.3 GeV and 
0.0021 < ycut < 0.0279 for data taken at 43.8 GeV. In figures 3 and 4 the hadron-level 
four-jet distributions with the fit results for the six energy points are shown together with 
the four-jet rate prediction with 0:3 being the current world average. The results of the 
fits are summarized in Table 2. The statistical uncertainties correspond to the uncer- 
tainty from the minimization, while the systematic errors are determined as described 
in section 4. 

It is of interest to combine the measurements of as from the different centre-of-mass 
energy points in order to determine a single value of as at a common energy scale. The 
fit results for as are combined using the procedure of Ref. [30]. In brief the method is as 
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[GcV] 


tts(v^) 


stat. 


exp. 


hadr. 


scale 


xVd.o.f. 


14.00 


0.1536 


0.0032 


0.0064 


0.0028 


0.0074 


1.46/3 


22.00 


0.1407 


0.0028 


0.0034 


0.0021 


0.0024 


14.22/6 


34.60 


0.1346 


0.0007 


0.0019 


0.0031 


0.0011 


17.20/7 


35.00 


0.1391 


0.0006 


0.0017 


0.0033 


0.0012 


23.51/7 


38.30 


0.1355 


0.0021 


0.0042 


0.0038 


0.0020 


19.78/8 


43.80 


0.1289 


0.0012 


0.0011 


0.0038 


0.0019 


4.02/9 



Tabic 2: The value of as for each energy point and the statistical, experimental, hadroniza- 
tion and scale errors. The last column corresponds to the x^/d.o.f. value of the fit. 

follows. The set of as measurements to be combined are first evolved to a common scale, 
Q = Mzo , assuming the validity of QCD and Mzo being the mass of the Z° vector boson. 
The measurements are then combined in a weighted mean, to minimize the between 
the combined value and the measurements. If the measured values evolved to Q = Mzo 
are denoted a;s,j, with covariance matrix V, the combined value, as{Mzo), is given by 

as{Mzo) ^^Wias,i where = — ^— — — ^, (6) 

where i and j denote the six individual results. The difficulty resides in making a reli- 
able estimate of V in the presence of dominant and highly correlated systematic errors. 
Small uncertainties in the estimation of these correlations can cause undesirable features 
such as negative weights. For this reason only experimental systematic errors assumed 
to be partially correlated between measurements contribute to the off-diagonal elements 
of the covariance matrix: V'- = mm[a'^^p ^, a^^^ j) . All error contributions (statistical, 
experimental, hadronization and scale uncertainty) contribute to the diagonal elements 
only. The hadronization and scale uncertainties are computed by combining the as val- 
ues obtained with the alternative hadronization models, and from the upper and lower 
theoretical errors, using the weights derived from the covariance matrix V. 
The fit result from the 14 GeV data has large experimental and theoretical uncertainties. 
We therefore choose to not include this result in the combination. The combination using 
all results for ^/s > 22 GeV is 

as{Mzo) - 0.1159 ± 0.0004(stat.) ± 0.0012(exp.) ± 0.0024(had.) ± 0.0007(theo.) , (7) 

consistent with the world average value of as{Mzo) — 0.1182 ± 0.0027 [31]. The weights 
Wi as described in Eq. 6 are 0.15 for 22 GeV, 0.29 for 34.6 GeV, 0.29 for 35 GeV, 0.06 for 
38.3 GeV and 0.21 for 44 GeV. The results at each energy point are shown in figure 5 and 
compared with the predicted running of as based QCD and on the world average value of 
as (Mzo). For clarity the values at 34.6 and 35 GeV have been combined at their luminosity 
weighted average energy ^/s — 34.8 GeV using the combination procedure described 
above. The combined value is as(34.8GeV) = 0.1368 ± 0.0005(stat.) ± 0.0017(exp.) ± 
0.0032(had.) ± 0.0011(theo.). The results of ALEPH [10], DELPHI [11] and OPAL [12] 
are shown as well. 
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5.3 Renormalization Scale Dependence of as 

For the fits in Section 5.2 the renormahzation scale factor is set to the natural choice 
X/^ — 1, where the energy dependence of Eq. 3 comes only from the running of the 
coupling as. However, different schemes for the determination of the renormalization 
scale factor are proposed. As a cross check of our default result we investigate in this 
section two more choices of the renormalization scale factor. 

In the optimized renormalization scheme [32] the minimization is performed for all energy 
points with as and treated as free parameters. The systematic uncertainties are 
determined using the optimized renormalization scale determined in the default fit. The 
value of as, the optimized scale and their correlation are summarized in Table 3. The 



A/i [GeV] 


as(v^) 


stat. 


exp. 


hadr. 




Corr. 


xVd.o.f. 


14.00 


0.1572 


0.0033 


0.0065 


0.0031 


1.49 ±0.50 


0.66 


0.01/2 


22.00 


0.1393 


0.0027 


0.0033 


0.0019 


0.55 ±0.31 


0.28 


13.38/5 


34.60 


0.1357 


0.0007 


0.0019 


0.0030 


0.44 ±0.07 


-0.68 


4.09/6 


35.00 


0.1401 


0.0006 


0.0017 


0.0031 


0.46 ±0.06 


-0.67 


6.04/6 


38.30 


0.1416 


0.0024 


0.0049 


0.0039 


0.33 ±0.05 


-0.68 


10.98/7 


43.80 


0.1291 


0.0012 


0.0011 


0.0037 


0.87 ±0.33 


-0.44 


3.90/8 



Table 3: The value of as and the statistical, experimental, hadronization, renormalization 
scale factor a;°P*, the correlation between as and the renormalization scale factor and the 
X^/d.o.f. value of the fit for all energy points with the minimization performed using the 
optimized renormalization scheme. 

variation of x^/d.o.f. as a function of the scale is shown in Figure 6. The combination 
of all energy points at ^/s > 22 GeV using the method described in Section 5.2 ^ yields 

as(Mzo) = 0.1161 ± 0.0005(stat.) ± 0.0012(exp.) ± 0.0022(had.). (8) 

The values for for centre-of-mass energies of 14. 22, 34.6, 35 and 43.8 GeV are within 
their statistical uncertainties covered by the systematic variation (0.5 < < 2.0) of the 
default fit. 

The second choice for the determination of the renormalization scale factor followed ap- 
proximately the approach of "minimal sensitivity" (PMS) suggested by [33]. The renor- 
malization scale factor x^™ is specified by the point with as having the least sensitivity 
to the renormalization scale factor x^. The variation of as a function of x^ is shown 
in Figure 6. The renormalization scale factor x™™ is determined by a fit to the as varia- 
tion with respect to x^ around the minimum. The determination of as is then repeated 
with x^ set to x™"^. At 14 GeV the variation of as with respect to the renormalization 
scale factor has no minimum and therefore no fit is performed. The renormalization scale 
factor a;™™ as well as the result of the fit are summarized in Table 4. The systematic 
uncertainties are determined using a;™™ determined in the default fit. The combination of 

^In this case the statistical, experimental and hadronization uncertainty only contribute to the diagonal 
elements of the covariance matrix V. 
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stat. 
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22.00 


0.1391 


0.0027 


0.0034 


0.0018 


0.42 


13.61/6 


34.60 


0.1345 


0.0007 


0.0019 


0.0031 


0.92 


15.38/7 


35.00 


0.1391 


0.0006 


0.0017 


0.0033 


0.92 


20.78/7 


38.30 


0.1354 


0.0021 


0.0042 


0.0038 


1.15 


20.80/8 


43.80 


0.1288 


0.0012 


0.0011 


0.0038 


1.22 


4.50/9 



Table 4: The value of as and the statistical, experimental, hadronization, the renormaliza- 
tion scale scale factor x^^^ and the x^/d.o.f. value of the fit for the energy points between 
22 and 43.8 GeV with as having the least sensitivity to the renormalization scale x^. 

all energy points above ^/s > 22 GeV using the method described in Section 5.2 ^ yields 

as{Mzo) = 0.1155 ± 0.0004(stat.) ± 0.0012(exp.) ± 0.0022(had.), (9) 

consistent with Eq. 7. The local minimum of as as a function of x^, is very close the 
natural scale a;^ = 1 leading to a fitted value of as similar to the default value Eq. 7. 
The choice of the renormalization scale factor x™"^ and a;°P* returns a value of as which is 
within the variation of the systematic uncertainty due to missing higher order terms. We 
therefore conclude that the evaluation of theoretical uncertainties using the commonly 
used standard method of setting x^ to 0.5 and 2. safely includes alternative methods of 
the choice of the renormalization scale factor. 



6 Summary 

In this paper we present measurements of the strong coupling from the four-jet rate at 
centre-of-mass energies between 14 and 44 GeV using data of the JADE experiment. The 
predictions of the PYTHIA, HERWIG and ARIADNE Monte Carlo models tuned by 
OPAL to LEP 1 data are found to be in agreement with the measured distributions. 
From a fit of QCD NLO predictions combined with resummed NLLA calculations with 
X/^—l to the four-jet rate corrected for experimental and hadronization effects we have 
determined the strong coupling as. In addition we investigated two more choices for the 
determination of the renormalization scale and found the results to be consistent. The 
value of as(Mzo) is determined to be as(^z") = 0.1159±0.0028 (total error). The natural 
choice of the renormalization scale a;^ = 1 is close to the renormalization scale factor with 
as having the least sensitivity to the renormahzation scale factor, ^ = 0. Therefore the 
theoretical uncertainty determined by setting the scale x^ to a;^=0.5 and x^=2.0 yields 
smaller theoretical uncertainties than e.g. for fits to event shape observables [16]. This is 
also true for measurements of as performed at LEP energies [12]. 
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Table 5: Hadron-level value of the four-jet fraction using the Durham algorithm at 14, 
22 and 34.6 GeV. In all cases the first quoted error indicates the statistical error while 
the second quoted error corresponds to the total experimental uncertainty. Uncertainties 
consistent with zero indicate that the corresponding value is smaller than the precision 
shown in the table. 
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Table 6: Hadron-level value of the four-jet fraction using the Durham algorithm at 35, 
38.3 and 43.8 GeV. In all cases the first quoted error indicates the statistical error while 

the second quoted error corresponds to the total experimental uncertainty. Uncertainties 
consistent with zero indicate that the corresponding value is smaller than the precision 
shown in the table. 
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Figure 1: The figures show the four- jet rate distribution corrected for detector effects 
at hadron-level as a function of the ycut resolution parameter obtained with the Durham 
algorithm. The four-jet rate at four average centre-of-mass energies are shown for ^/s = 14 
to 35 GeV in comparison with predictions based on PYTHIA, HERWIG and ARIADNE 
Monte Carlo events. The expected bb background is subtracted from the data. The error 
bars show the statistical (inner part) and the experimental and statistical uncertainties 
added in quadrature. The panel in each upper right corner shows the differences between 
data and Monte Carlo, divided by the quadratic sum of the statistical and experimental 
error. At points with no data events, the difference is set to zero. 
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Figure 3: The plots show the four-jet rate distributions at the hadron-level for ^/s = 
14 GeV and 22 GeV. The error bars show the statistical (inner part) and the experimen- 
tal and statistical uncertainties added in quadrature. The solid curve shows the theory 
prediction after minimization within the fit range indicated. The dash-dotted lines 
shows the error band of the four-jet rate prediction with as(Mzo) being the current world 
average value and its error [31]. 
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Figure 5: The values for as at the various energy points. The errors show the statistical 
(inner part) and the total error. The full and dash-dotted lines indicate the current world 
average value of as(-^zo) [31] with error. The results at y/s — 34.6 and 35 GeV have been 
combined for clarity. The results from the LEP experiments ALEPH [10], DELPHI [11] 
and OPAL [12] are shown as well. 
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Figure 6: The result of as(Mzo) and the x^/ d.o.f. of the fit to the four-jet rate as a function 
of the renormahzation scale for ^/s = 14 GeV to 43.8 GeV. The arrows indicate the 
variation of the renormahzation scale factor used for the determination of the systematic 
uncertainties. 
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